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RANDOM AND SYSTEMATIC MEASUREMENT ERRORS 


IN ACOUSTIC IMPEDANCE AS DETERMINED 
BY THE TRANSMISSION LINE METHOD 

Tony L. Parrott and C. D. Smith* 

Langley Research Center 

SUMMARY 

The effect of random and systematic errors associated with the measurement of 
normal incidence acoustic impedance in a zero -mean -flow environment has been inves- 
tigated by the transmission line method. The investigation was both experimental and 
analytical. The influence of random measurement errors in the reflection coefficients 
and pressure minima positions was investigated by computing fractional standard devia- 
tions of the normalized impedance. Both the standard techniques of random process 
theory and a simplified technique were used. 

Over a wavelength range of 68 to 10 cm (frequency range, 0.5 to 3.5 kHz), random 
measurement errors in the reflection coefficients and pressure minima positions could 
be described adequately by normal probability distributions with standard deviations of 
0.001 and 0.0098 cm, respectively. The assumption of normal probability distributions 
for the measurement errors simplified the application of standard techniques in the cal- 
culation of standard deviations in the impedance components. An error propagation 
technique based on the observed concentration of the probability density functions was 
found to give essentially the same results but with a computation time of about 1 percent 
of that required for the standard technique. 

The results suggest that careful experimental design reduces the effect of random 
measurement errors to insignificant levels for moderate ranges of test specimen imped- 
ance component magnitudes. Most of the observed random scatter can be attributed to 
lack of control by the mounting arrangement over mechanical boundary conditions of 
the test specimen. 


INTRODUCTION 

Measurements of normal incidence specific acoustic impedance continue to be 
important in determining the acoustic properties of materials considered for use in 
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noise control. For materials that can be classed as locally reacting, the normal incidence 
impedance can specify the boundary condition in any system amenable to analysis. For 
homogeneous and isotropic materials characterized by extended reaction, two sets of 
normal incidence impedance measurements are sufficient to determine the material prop- 
agation constant and characteristic impedance. (See ref. 1.) Recent investigations 
dealing with the optimum design of noise -attenuating lined ducts (ref. 2) have raised ques- 
tions about the degree of precision and accuracy needed to specify duct liner impedances. 
Reference 3 gives in situ impedance measurements on a duct liner at 20 axial locations 
2.54 cm apart. These measurements indicated as much as a 30-percent random scatter 
about the mean trend. A limited error analysis suggested that the observed scatter was 
dominated by material inhomogeneity rather than by impedance measurement errors; 
however, a detailed discussion of the effects of random measurement errors was not 
given. 

This investigation isolated the effects of measurement errors in the determination 
of acoustic impedance apart from the effects of material inhomogeneity. The investiga- 
tion was limited to thin face-sheet materials studied in a zero-mean-flow environment 
by means of the transmission line method. Three different test specimens (designated 
A, B, and C) were used in the investigation. Specimen A was chosen because its imped- 
ance could be predicted theoretically. Therefore, a comparison of the theoretical and 
measured impedance for specimen A allowed the combined effect of all systematic errors 
to be evaluated. Test specimens B and C were fabricated from fiber metal. Repetitive 
measurements on all the specimens over the frequency range of 0.5 to 3.5 kHz in 
0.1-kHz increments provided sufficient data to establish meaningful statistics for ran- 
dom variations in the measured reflection coefficients and pressure minima positions. 

Two methods were used to calculate the standard deviations in the deduced test 
specimen impedance from the statistics of the measurements. In one method the stand- 
ard techniques of stochastic process theory (hereafter called the exact method) were 
used to calculate the variances of the impedance components in terms of the functional 
relation between the impedance and the directly measured quantities. The second method 
(hereafter called the approximate method) resulted from a series expansion of the func- 
tional relation connecting impedance and the directly measured quantities. The resulting 
approximate formulas involved the first-order derivatives of the impedance components 
with respect to the measured quantities and their variances. 

SYMBOLS 

A = ^ cm®, or ^ cm^ 

i=l i=l 


2 


amplitude of acoustic wave incident on test specimen, dynes-cm 


-2 


B 


^0 

'p 

f 

SPO 

fx(^) 

JoO 

■>2<) 

j 

k 


amplitude of acoustic wave transmitted through test specimen, dynes-cm 
N 

= ^ n., cmO, or ^ x., 


-2 


cm 


i=l i=l 

amplitude of acoustic wave reflected from test specimen, dynes-cm 

-2 

amplitude of acoustic wave reflected from cavity back, dynes-cm 

-1 


-2 


sound speed at temperature T, cm- sec 


-1 


sound speed at reference temperature Tq, cm-sec 
specific heat of air at constant pressure, cal-g"^-°C ^ 


error propagation functions 


frequency, Hz 


probability density function for random variable 


probability density function for random variable Xj 
probability density function for random variable 


probability density function for random variable Rq 


probability density function for random variable X 


zeroth-order Bessel function 


second-order Bessel function 


imaginary unit 


wave number, 2tt/X, cm 


-1 


3 



I channel length, cm 

M arbitrary integer 

N total number of pressure minima 

n complex polytropic gas constant 

n^ ith pressure minimum 

P porosity 

Pmax maximum acoustic pressure, dynes-cm“^ 

p . minimum acoustic pressure, dynes-cm"^ 

^min ^ 

-2 -1 

R specific acoustic resistance, g-cm -sec 

reflection coefficient at ith pressure minimum 
Rq reflection coefficient at test specimen surface 

r radial coordinate, cm 

r^ channel radius (test specimen A), cm 

rp probe radius, cm 

rQ reciprocal of reflection coefficient 

standing -wave ratio associated with ith pressure minimum, dB 
Sq standing -wave ratio at specimen surface, dB 


s 


shear wave number. 
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T 

To 

X 


Xv 


Z 


a 

/3 

r 

y 

■n 

6 

'^2 

X 


workii^ temperature, °C 


reference temperature, 

-2 -1 

specific acoustic reactance, g-cm -sec 


tube axial coordinate, cm 


backing cavity depth, cm 


position of ith pressure minimum relative to test specimen surface, cm 

-2 -1 

specific acoustic impedance at test specimen surface, g-cm -sec 

-2 -1 

specific impedance of backing cavity, g-cm -sec 

-2 -1 

complex characteristic impedance of test specimen channel, g-cm -sec 


tube-wall absorption coefficient, cm 


-1 


square root of Prandtl number. 




complex propagation constant for single channel of test specimen A, cm 


ratio of specific heats for air 


-1 


specific acoustic impedance ratio, Z/pc 
nondimensional radial coordinate, 
acoustic resistance ratio, R/pc 


thermal conductivity, cal-cm’’^-sec"^-^C”^ 


parameter in expression for f 
acoustic wavelength, cm 
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Ml 


M2 




viscosity of air, g-cm"^-sec 


parameter in expression for f^ 

^1 


parameter in expression for 

^1 



p density of air, g-cm 

a standard deviation of parameter indicated by subscript (for example, is 

standard deviation of R^) ^ 

2 2 
(j variance of parameter indicated by subscript (for example, is vari- 
ance of Rp ^ 

2 2 
(7 variance of distance to ith pressure minimum, cm 

^i 

X acoustic reactance ratio, X/pc 

O- ^ correlation coefficient between best estimates of x^ and \ 

(V angular frequency, rad- sec ^ 

Subscripts: 


i pressure minima positions counted from test specimen surface 

m mean value of random variable 

0 value of quantity at test specimen surface or fixed value of radial coordinate 

Abbreviations: 

cgs centimeter -gram-second system 

dB decibels referenced to 0.0002 dynes/cm^ where 1 dyne = 0.00001 N 

h.o.t. higher order terms 
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I.D. inside diameter 

O.D. outside diameter 

SPL sound pressure level, dB 

Special notations; 

erf( ) error function 

" best estimates obtained from linear regression curve 

normalization by wavelength X 

SYSTEMATIC AND RANDOM ERRORS IN ACOUSTIC 
IMPEDANCE MEASUREMENTS 

The specific acoustic impedance at the surface of a homogeneous material is under- 
stood to mean the ratio of the acoustic pressure to particle velocity suitably averaged 
over a representative area of the material. Since the direct measurement of acoustic 
particle velocity is not practical, acoustic impedance must be measured by indirect 
means. The acoustic phenomenon central to the transmission line method of measuring 
impedance is a one -dimensional standing-wave field. This standing-wave field is gen- 
erated in a rigid wall tube by a source at one end and a sample of the test material at the 
opposite end. For an inviscid acoustic medium, the resulting standing -wave field wovdd 
be characterized by equal pressure maxima and equal pressure minima located at half- 
wavelength intervals along the normal to the material surface. The acoustic impedance 
can be deduced from measurements of the ratio of the pressure maxima to the pressure 
minima, the distance from the material surface to the first minimum, and the acoustic 
wavelength. 

Influences of the measurement apparatus and associated instrumentation on the 
acoustic field and test specimen behavior cause systematic errors in the deduced imped- 
ance. Also, there are random errors induced in the directly measured quantities due to 
uncontrollable environmental and/or operational factors. The effect and control of sys- 
tematic errors has been discussed in references 4, 5, and 6. This paper, therefore, 
concentrates on the effect of random measurement errors on the deduced impedance; 
however, systematic error sources relevant to this work are discussed briefly. 

Systematic Errors 

Systematic errors associated with the transmission line method for determining 
impedance can be conveniently divided into the three following categories: 
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(a) Errors intrinsic to impedance tube operation 

(b) Errors intrinsic to test specimen mounting 

(c) Errors intrinsic to the measurement instrumentation 

References 4 and 5 give detailed accounts of systematic errors intrinsic to the imped- 
ance tube operation. Therefore, only the results relevant to this investigation are dis- 
cussed here. A discussion of systematic errors intrinsic to the test specimen mounting 
and measurement instrumentation is contained in the section entitled "Experimental 
Setup and Test Procedure.” 

Effects of tube attenuation.- Figure 1 shows an impedance tube setup for deter- 
mining the impedance across thin face-sheet materials. A variable cavity depth xjj 
provides a known impedance for the wave A 2 transmitted through the test specimen. 
The interaction of incident and reflected waves Aj and B^, respectively, produces a 
standing -wave field. When measured by an axial pressure probe, a typical standing- 
wave field appears as shown in the sketch of figure 1. In this sketch the coordinate sys- 
tem origin is taken at the face of the test specimen with the positive axis running toward 
the source. Successive standing-wave ratios ^i.e., Pmax/^min) minima 

positions are denoted by Sj and x^, respectively, where i = 1, 2, . . ., N. Note that 
Sq denotes the ratio p 
measurable . 




max/"^min 


at the specimen surface which is not directly 


Fluid viscosity and thermal conduction in the acoustic boundary layer near the tube 
wall attenuate the wave motion (refs. 4 and 5). As figure 1 illustrates, this attenuation 
creates an increase in the pressure minima (i.e., standing -wave ratios decrease) with 
distance from the test specimen. Therefore, the standing -wave ratios Sj^, which result 
mainly from the test specimen-cavity system, are contaminated by tube -wall attenuation. 
The resulting difference between the standing-wave ratios Sq and is usually less 
than 1 dB; however, the use of Sj to approximate Sq results in significant systematic 
errors for some wavelengths and test specimen impedances. 

The systematic error Sq - Sj can easily be corrected by conversion of the mea- 
sured Sj to equivalent reflection coefficients which approximate a straight line as 
illustrated in the lower sketch of figure 2. A linear regression curve through the mea- 
sured Ri then provides an estimate Rq of the reflection coefficient Rq. The formula 
for converting the measured S^ to equivalent Rj, where Sj are measured in dB, is 
given by 


Ri = 


10 


5^/20 _ ^ 


S./20 

10 +1 


( 1 ) 
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A secondary effect of tube -wall attenuation is a small shift of the pressure minima 
positions toward the specimen surface. In terms of acoustic wavelengths the expected 
shift in the first pressure minimum position according to reference 5 is 




( 2 ) 


At 1.0 kHz, typical values of the parameters in equation (2) are 

^ = 0.0017 
k 

Xj = 0.25 cm 

Rj = 0.52 (Si = 10 dB) 


Therefore, 


Axi a 0.0001 cm 

This result suggests that systematic shifts in pressure minima due to tube-wall attenua- 
tion are relatively insignificant except in those cases where Ri becomes very small. 

In such cases, the pressure minima are not well defined because of the very small 
standing -wave ratios. Consequently, random errors would remain dominant. 

E ffect of backing cavity .- The impedance at the test specimen surface, as deduced 
from standing -wave measurements, is determined by the interaction of the impedance 
across the test specimen and the cavity impedance as illustrated in figure 1. For thin 
test specimens (i.e., thickness small relative to the wavelength) the total impedance at 
X = 0 is simply the sum of the specimen impedance and the cavity impedance. To avoid 
large systematic errors in the deduced impedance across the test specimen, the cavity 
depth must be adjusted for zero reactance. The cavity impedance is given by 




cot kx^ 


Therefore, the appropriate cavity depth is given by 


(3) 
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(4) 


where f is the driving frequency, T is the cavity gas temperature, a is the tube- 
wall absorption coefficient, k is the wave number, Cq is a reference sound speed, 
and Tq is the corresponding reference temperature. 

Effect of finite probe size.- A systematic error in the measurement of pressure 
minima positions can arise because of the finite size of the axial pressure probe. From 
reciprocity considerations, the Rayleigh end correction of 0.85rp for an unflanged pipe 
radiating into free space should be on the order of the probe end correction. However, 
the geometry of the probe end and of the surrounding hardware boundaries deviates sig- 
nificantly from the geometry for the Rayleigh end correction. Therefore, an experi- 
mental determination of the end correction is necessary. The experimental procedure 
for accomplishing this correction is discussed in the section entitled "Experimental 
Setup and Test Procedure." 

Effect of structure -borne vibration.- Mechanical vibration of the impedance tube 
walls can cause serious systematic disturbances of the standing -wave field. In addition, 
the excitation of the pressure transducer by structure -borne vibrations can mask low 
sound pressure levels at pressure minima when the standing -wave ratios are large. 
Experimental design procedures to minimize these sources of systematic error are dis- 
cussed further in the section entitled "Experimental Setup and Test Procedure." 


Random Errors 

The formulas for calculating the specific acoustic impedance at the surface of a 
test specimen from standing -wave data are well known (refs. 4, 5, and 6). These form- 
ulas, in slightly modified form for application in this work, are given for the resistance 
and reactance ratios, respectively, as 


e = 



1 + 


^1 - R| 
,1 + Rr 


1 + tan'^ 277X, 


^ ' tan^ 27TX, 


(5a) 


10 


V 



and 


X 


1 - 


1 - Rr 


il + R 


•0/ J 


tan 27TX, 


1 + 


1 -Rc 

.1 + 


tan 2 ttxi 


(5b) 


where Xj = The quantities and X are best estimates of the first pressure 

minimum position and the wavelength, respectively. The term Rq is a best estimate 
of the reflection coefficient at the specimen surface. The specific acoustic impedance 
ratio ^ at the test specimen surface is given by 




(6» + jx) 


(6a) 


where Z is the specific acoustic impedance at the specimen surface. Since — « 1 

k 

for most measurement situations, it follows that 


^ » 0 + jx (6b) 

The effects of random measurement errors on Rq and Xj and therefore on the 

impedance components 6 and x were analyzed by first calculating the variances 

and a 2 in terms of the variances a*. ^ and a- The variances a?, 2 2 

X Rq Xj Rq Xj 

were then calculated from the estimated variances ^ and a respectively. The 

Ki X. 

2 2 

computations for cr^ and cr^ were accomplished using a standard technique from 
random process theory and a simplified approximate method, each of which is outlined. 

Exact method .- The following assumptions were used to calculate and ct ^ 

by the exact method: 

(1) The probability density functions f#, 

Kq 

tistically independent. 

(2) The density functions fj^ and f^^ ^xA are normally distributed. 


(Rq) and fc^(5j) 


are treated as sta- 


ll 



(3) If the respective means are suppressed, then 


‘xiW ='x. 


Therefore, from the random process theory given in reference 7 


l,m x-|^ P U,m R| 


« P i,m X;. p u,m ^ / v \ ^ 

a = L (^o) (^l) 

® -Jx. -Mai ^Rn,^-Mae, ^ ^ 

l,m Xh U,m Rrt 


Xi ^+Mai ™+Ma^ 

i,m x-i p U,m R 


2 p i,m Xj p u,m 2 /- \ /- \ - t 

fX - Xrv.1 Rn) 4 xJ dRn dx. 

X Jx, .. -Mai Jr. ...-Mai, i RqV O; x^V V 0 3 




(7b) 


where the mean values 0^^ and given by 


m 


= 1. 


Xi _+Map, Rn _+Maf 


i,m Xj p u,m -^0 o \ ^ 

Jr. 




(8a) 


Xi .,+Mai R. ^+Man 
1 , m X *1 p 0 , m Xv| 


^m 


p x,m X| p u,m n,Q . v /^ \ - 

= u " xfp (Ro) ^5 (^i) 

Jx, „-Ma?. Jr. „.-Ma- ^ RgV 0/ XjV ly 0 1 


-'x, ,^-MaJ ^R. -Ma„ 
l,m x^ u,m "^Q 


(8b) 


The integration limits in equations (7) and (8) were truncated symmetrically about the 

mean values R. ^ and x, ^ at integer multiples of aA and a^ , respectively, by 
u , m j. ^ni ^0 1 

trial and error. Values of cr^ and a$ were derived from estimates of ^ 

^ 2 ^0 ^1 “i 

and a 

H 

The variances a„ ^ and a ^ were estimated by compiling statistical frequency 

X| 

distributions of the deviations of measured values of Rj and x. from corresponding 
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best estimates obtained from their respective linear regression curves. A typical 
regression curve for four measured reflection coefficients (i.e., N = 4) is shown in the 
sketch of figure 2. The best estimate of the specimen reflection coefficient is the extrap- 
olated value Rq obtained from the regression curve at x = 0. Likewise Rj are the 
best estimates of Rj^ taken from the regression curve at the corresponding x^. The 
deviations (r^ - R.) are shown in figure 3 as a probability density histogram. Super- 
imposed on figure 3 is a normal probability density function adjusted so that the inte- 
grated areas are equal between the deviations -0.003 and 0.003. Therefore, an estimate 
of f (Rj) is given by 


f 




exp 



(9) 


where o-t, * 0.001. 

Ri 

An entirely similar procedure for estimating the density functions 
trated in figures 4 and 5. In this case the result is 




is illus- 




exp 




2a 2 
^i 


( 10 ) 


where a « 0.0096 cm. 

^i 

To evaluate the integrals in equations (7) and (8), values of 
needed to define the integration limits and to construct f|^ (^o) 



and 


and ai 



are 

which 


appear in the integrands. The random variable Rq is linearly dependent on the mea- 
sured random variables R^ through the linear regression curve shown in figure 2. 
Therefore, from the equation for the linear regression curve and the definition of vari- 
ance, it follows that 



( 11 ) 
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where 


N 

= 


i=l 


N 
i=l 


B= > X. 


and 


N ^ 2 


The integer N denotes the total number of and corresponding measured for a 
given wavelength. 

The procedure for calculating ^ from ^ is straightforward but lengthy. 


^i 


The procedure is outlined and the results are stated. The first step is to calculate 
a- 2 and in terms of The same procedure used to obtain equation (11) 


*1 

gives 


^i 


a-. 2 = ^ 2 /A-2B-H N' 
^i ^i V NA - b2 


(12a) 


2 _ „ 2 2N 
^ ^ ^i 


(12b) 


where 


N 

i=l 


B 


N 


i=l 


and 


N 22 
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The behavior of the ratios a- ^ and acVa ^ is shown by the respective curves 

xi/ X. X. 

2 2 

of figures 6 and 7. These curves show that both and decrease with 


2 2 
increasing N; however, decreases at a much greater rate than does . 


Although the density functions fx.(xi) are assumed normal, the density 


function 


IS non- 


fi (xi) will not be normal since the random variable transformation Xj^ = i 

linear. The technique for deriving from f^^^Xj^^ and f^(X) is given in ref- 


erence 7. The result is 




Xl,A 


exp( — - Kgj 
fig 


expf-iq 

\ ^2/ ^ \JJ^2 


(13a) 


where 


V = 


2 (1 - nf 

XfjXI Xj^ A 


^\\l f - - 

^x^,X^m “ ^l,m ~ ^m 


(13b) 


Mo = 


211 -"IJV 


\ 2 /„ V 2 

CT- \ /a- \ 

Xj\ / x,\ / 2 

1 x,| - j I + I I [ 1 “ ^0 X 

1 Xj,X J \^x \ ^1’''- 


(13c) 


'^2 = 




a- /(Tr- ' ^ 


IX, „ - n- cX 


l,m X, ,X m ac 


+ X 


2 

m \ at 


1 - n: 


Xi,X 


(13d) 


The quantities x, and X^ are the true mean values of the random variables 
^ l,m m 

Xj^ and X. In the data gathering process only one set of x^ values was obtained for 
each wavelength on a given test run. Therefore, Xj^ ^ and X^ were estimated 
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is given by 


by Xj and X. An expression for the correlation coefficient ^ 

N - B 


a 


\/nA - b2 \/an - 2BN + 


( 14 ) 


The behavior of <• as a function of N is shown in figure 8. As would be expected, 

Xj^,X 

^ decreases with increasing N. 

2 2 

Equations (5) and (7) to (15) permit exact calculations of and from 

experimentally determined input parameters. A computer program was prepared to 
perform these calculations for several representative values of the various parameters. 
Because a large amount of central processing xmit (CPU) time (typically, 100 sec) was 
required to perform a trend analysis, a simplified approximate analysis was developed. 
The results of the exact and approximate analysis are compared in this paper. However, 
before discussing these comparisons, the approximate analysis procedure is outlined. 


Approximate an alysis.- The approximate analysis follows the procedure given in 
reference 7 which makes use of the assumptions 


(a) The density functions f|^ 

the vicinity of Rn ^ and x, 

•' u,m l,m 


and f; 


ii(%) 


are smooth and concentrated in 


of R, 


(b) The impedance components 0 ^Rq,Xj^ and x(Ro>^l) smooth in the vicinity 


0,m 


and X 


l,m 


and do not take on large values out of the neighborhood of R, 


0,m 


and X 


l,m* 


(c) The random variables Rq and Xj^ are uncorrelated. 

The first terms of the Taylor series expansion of the resistance ratio about the mean 


values R, 


0,m’''l, 


m 


give 


«(Ml) » 


R, 


0,m’^l,: 


(«0 - ^, m ) * 


m 


R, 


0,m’^l, 


1 l,ni 


,m 


+ h.o.t. 


Therefore, 
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, /o i \ 


- H^0,m’^l,mj 



2 




2 


+ 


30 ' 


L 3x 




m 




2 


+ 



Substituting equation (15) into the integral of equation (7) and usii^ assumptions (a), (b), 
and (c) yield 


a 





(16a) 


In the same manner, for the reactance ratio 


CT 


X 


2 



(16b) 


where the derivatives in equations (16) are evaluated at the appropriate mean values 

Flo and x-. . Equations (16) are general forms relating the variance of a dependent 

vi j in ijiii 

variable to the variances of any number of independent variables. Therefore, the vari- 

p 

ance crt can be written 
^1 



It is convenient to change independent variables in equation (16) as follows; 


^0 ■ 


1 + Rq 


(18a) 
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(18b) 



= 2irx^ 


Therefore, 


^0 


2 







where 


cr 


2 

0 







(19a) 


(19b) 


(20a) 


(20b) 


If the derivatives in equations (19a) and (19b) are calculated using equations (5), the fol- 
lowing matrix equation can be written: 




r F 2 
^2 



( 

CSI 

o 

b 

\ 




F 2 
^2 




where 


2 ^ 
sec kxj^ tan kx^ 

^1 + rQ^ tan^ 


(21a) 


(21b) 
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(21c) 


9^ / >k2 2-^\ 

sec'^ kxj (1 - Tq tan kxj^j 

^ / •‘O 

f 1 + Tq 2 tan^ kxj^j 

2 2 

Equations (18) and (21) are sufficient to calculate and cr^ in terms of the experi- 

mentally determined parameters. The behavior of the error propagation functions Fj 
and Fg, as Rq and range through values likely to be encountered in typical imped- 

ance measurements, should provide insight as to which regions in the measurement 
space may be critical in regard to random measurement errors. Figures 9 and 10 show 
plots of Fj^ and Fg, respectively, as functions of with rQ as a parameter. Note 
that rQ can be physically interpreted as the reciprocal of the standing-wave ratio at the 
specimen surface, and kxj has been written as 27 txj^ to emphasize the fact that the 
distance to the first pressure minimum is measured in wavelengths. The parameter rg 
extends from 0.2 to 0.5, a range which corresponds to a standing -wave ratio range from 
6.0 to 14.0 dB or a normal incidence absorption coefficient rar^e from 0.56 to 0.89. The 
standing -wave ratios in dB corresponding to the values of Vq are also shown in fig- 
ure 9. Since the functions Fj and Fg are periodic with a period of n radians, their 
complete behavior is described over the Xj interval 0.00 to 0.50. Also, since the 
squares of Fj and Fg enter equation (21a), the behavior of the variances and 

is symmetrical about x^^ = 0.25. The plots of figures 9 and 10 clearly indicate 
that random measurement errors tend to play an increasingly important role as Xj 
approaches 0.25 and as the standing-wave ratio increases beyond about 6.0 dB. This 
result corresponds to what would be expected from physical considerations because the 
acoustic particle velocity tends toward 0 as the impedance at the specimen surface 
increases indefinitely. 

Comparison of exact and approximate analysis. - It remains to be seen whether the 

approximate analysis yields sufficiently accurate values of a„ and ct given the 

^ X 

restrictive assumptions underlying the analysis. Because computation time required in 
the approximate analysis is almost a hundredfold less than the time required in the exact 
analysis, a comparison of the results of the two analyses for some t3qjical input param- 
eters is worthwhile. The results of this comparison are shown in the curves presented 
in figures 11 to 16. In these figures, the solid curves represent the results from the 
approximate analysis and the dashed curves represent results from the exact analysis. 
Each set (parts (a) and (b)) of curves in figures 11 to 14 present computed values of the 
fractional standard deviations in nondimensional resistance ^■nd in nondimensional 

reactance Plotted against the reflection coefficient Rq with the nondimensional 

distance to the first pressure minimum Xj as a parameter. Realistic assumed values 
for the standard deviations in the reflection coefficients, the pressure minima positions. 
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and the acoustic wavelengths were chosen for the computations. Preliminary measure- 
ments suggested that the standard deviations in the reflection coefficients could be taken 
as approximately 0.002 and the standard deviations in the pressure minima positions 
could be taken to be approximately 0.02 cm. The four sets of computations were made 
for wavelengths corresponding to frequencies of 500, 1000, 2000, and 3000 Hz. 

Figures 11(a), 12(a), 13(a), and 14(a) show that for Xj^ = 0.1 (i.e., provided the 
pressure minimum position is not in the immediate neighborhood of the quarter-wave 
point) the fractional standard deviations in the resistance Ogy'e increase relatively 
slowly as Rq increases up to a value of about 0.7. Beyond Rq = 0.7, however, 

increases rapidly and can exceed 6 percent for values of Rq greater than about 
0.975, a value which corresponds to a normal absorption coefficient of about 5 percent. 

As the pressure minimum position approaches the neighborhood of Xj = 0.25, 
begins to increase at smaller values of Rq. Furthermore, the increase becomes more 
pronounced with increasing frequency as indicated by the curves for Xj = 0.200, 0.225, 
and 0.245 (note scale change in fig. 13(a)). For the values of o and chosen, 

Xi Kj 

the exact calculations (dashed curves) agree well with the approximate calculations. The 
discrepancy that does exist could result from the assumption in the approximate calcula- 
tion that the probability density function is concentrated at the mean values. 

Figures 11(b), 12(b), 13(b), and 14(b) present the same parametric study of the 
fractional standard deviation of the reactance ratio as was done in figures 11(a), 12(a), 
13(a), and 14(a) for the resistance ratio. These figures show that more 

sensitive to values of x^ than was Also, the behavior of ^^X with increasing 

Rq is in general more complicated than was the behavior of Again, because of 

the periodicity in the trigonometric functions involved in Fj^ and F 2 , the curves for 
Xj >0.25 repeat the behavior of those for x^ <0.25 and, therefore, were not plotted. 
Although the exact and approximate calculations agree well for the cases presented, the 
discrepancy between the two methods of analysis is generally greater for a N than 
for 

In addition to the comparisons of the approximate and exact analysis presented in 
figures 11 to 14, a comparison of the two methods of analysis for increasing values of 
a-p and a also yields interesting results. Figures 15 and 16 illustrate the results 

JrCi Xj 

of this comparison for a frequency of 3.0 kHz. Figures 15(a) and (b) show the behavior 
of the fractional standard deviations in resistance and in reactance ratios, respectively, 
as a function of the standard deviation in reflection coefficients with quantities held con- 
stant as denoted in the figure key. As in the previous comparisons, the approximate 
calculations are shown by the solid lines, and the exact calculations are shown by the 
dashed lines with Xj as a parameter. Clearly, the approximate calculation begins to 
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break down for values of 


CTd exceeding approximately 0.006 for both the fractional 
”i 


standard deviation in the resistance and in the reactance ratio. Figures 16(a) and (b) 
show the behavior of the fractional standard deviations in resistance and in reactance 
ratios, respectively, as a function of the standard deviation in the pressure minima posi- 
tions, with quantities held constant as denoted in the figure key and with Xj^ as a param- 
eter. In this case, the approximate calculation begins to break down for values of ct 
exceeding approximately 0.045 cm. 


Figures 11 to 16 are believed to contain fairly representative comparisons of the 
approximate and exact methods for calculatir^ the standard deviations of the impedance 
components when given the standard deviations in the reflection factors and in the pres- 
sm*e minima positions. In summary, the following observations are noted; 


(1) Fractional standard deviations in resistance ratios tend to increase with 
increasing reflection coefficient. 

(2) Fractional standard deviations in reactance ratios tend to develop complicated 
dependencies on the reflection coefficient as the frequency increases and as the pressure 
minima positions approach the quarter -wave length point. 


(3) Assumption (a) in the approximate analysis procedure appears to break down for 

values of and a exceeding about 0.006 and 0.045 cm, respectively. 

«i Xj 

(4) For values of and a below 0.006 and 0.045 cm, respectively, there is 

^i ^i 

generally good agreement between the approximate and exact procedures; however, the 
discrepancy does increase with frequency and as Xj approaches the quarter-wavelength 
point. Also, the discrepancy is greater between the exact and approximate methods for 



These results indicate that for the transmission line method of impedance measure- 
ments, the simple approximate equations (20) and (21) are suitable for providing conser- 
vative estimates for the variances in the impedance components in terms of the variances 
of the measured quantities Or and provided they do not exceed 0.006 and 0.045 cm, 

XVi 

9 9 

respectively. For a given test specimen, ^.nd are dependent on the quality 

of the particular experimental test setup and instrumentation, variability of environmental 
parameters such as temperature and ambient acoustic noise, and mechanical vibration. 


Methods for estimating vary from mere guesses to 

elaborate statistical procedures. In this investigation these estimates were obtained by 
means of a simple compilation of statistical data associated with the linear regression 
curves for obtaining the estimates R^, Xj, and X as described in the section entitled 
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"Systematic Errors," The next section describes the data collection and processing in 
more detail. 


DATA COLLECTION AND PROCESSING PROCEDURE 


The random measurement errors associated with reflection factors R^ and pres- 
sure minima positions can be quantified by constructir^ parametric probability den- 


sllies t. 


) 


‘XjV"l’*'i m’ 'm/ (^i’^i m’^m) ®^ch of a series of test specimens 

selected to produce a range of reflection coefficients between 0 and 1,0, Such a procedure 
wovild require several hundred repetitive measurements for each specimen and for each 
test waveler^th. For example, 5 different reflection factors and 5 different wavelengths 
would imply 25 probability density functions. Experience with standing -wave measure- 


ments suggests that the dependency between (^i’^i 


m’^m; 


■i’"i,m’^m 


.) 


IS 


probably very weak except in extreme and, from a practical standpoint, uninteresting 
situations where the test specimen is highly absorptive (i.e., R^ — O) or highly reflec- 
tive (i.e,, R| — l). For the highly absorptive case, the standing-wave ratios are small 
(typically less than 1 dB); therefore, the random errors in pressure minima positions are 
large because the pressure minima positions are not well defined. On the other hand, 
when the reflection factors are large, the pressure minima positions are well defined, 
but the standing -wave ratios are large and tend to be unstable (i.e., subject to environ- 
mental disturbances). In addition, the acoustic wavelength also affects the statistical 
dependency between the random variables x^ and R.. 

The inordinate number of measurements required to rigorously construct the prob- 
ability densities ^x. (^i’^i m’^m) ^R (^i’^i m’^m) can be dramatically reduced 

with the help of the following assumptions: 


(1) Random measurement errors in reflection factors and pressure minima posi- 
tions can be treated as statistically independent if the test specimen is not too highly 
absorptive or reflective. 


(2) Deviations of individual measurements from their respective parent population 
means can be estimated by the deviations from their respective best estimates as 
obtained from linear regression curves (figs. 2 and 4) for a given impedance 
measimement. 


Assumption (1) appears plausible for the range of standing -wave ratios and wave- 
lengths typically encountered in impedance measurements of most interest. The statis- 
tical bias resulting from assumption (1) will tend to cause overestimates of 2 and 
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2 

or whereas the statistical bias resulting from assumption (2) will tend to cause under- 

Xi 

estimates of a-r, ^ and a 2, 

Ri X. 


Data Collection 

To generate the deviations - R. and x^ - x^, at least three standing-wave 
ratios and corresponding pressure minima positions x^ were measured for each 
test frequency. For the three test specimens designated A, B, and C, the total number 
of deviations measured exceeded 3500. A set of deviation measurements for a given 
specimen was collected over a frequency range of 0.5 to 3.5 kHz at typically 0.1-kHz 
intervals. Standing -wave ratios could be read in dB to three significant figures (i.e., 
±00.1 dB) and pressure minima positions in cm to four significant figures (i.e., 

±00.01 cm). Approximately 4 hr were required for collecting a data set. During that 
time the ambient temperature varied by approximately 2® C. From three to five data 
sets were collected for each specimen. As much as 6 months elapsed between some 
data sets. Except for test specimen A, the test specimens were removed and reinstalled 
in the test fixture between test runs. 


Data Processing 

A computer program was written to process the standing-wave ratio and pressure 
minima data. The computer program involved the following calculation steps: 

(1) Conversion of to R. 

(2) Calculation of linear regression curve coefficients for R^ and x^ 

(3) Calculation of the best estimates R^^, x^ X, x^^, and Rq 

(4) Tabulation of x^ - x^ and R^^ - R^ into selected class intervals 

(5) Calculation of 6 and x from Rq and x^ 

(6) Calculation of fractional standard deviations and from o’ 

and a ^ using the approximate analysis results given by equations (18) 
i 

and (21) 

2 2 

The estimates of a„ ^ and a ^ in step (6) undergo improvement as more data are 

Ki x^ 

collected and tabulated according to step (4). After sufficient data are collected, the 
estimates ^ and as shown, respectively, in figures 3 and 5, change little pro- 

vided that operational and environmental parameters are held constant. 
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EXPERIMENTAL SETUP AND TEST PROCEDURE 


Two equally important objectives were considered when the experimental setup 
and test procedure were designed. One objective was to expedite the accumulation of 
standing -wave impedance data by readii^ standing -wave ratios and pressure minima 
positions directly from the monitoring instrumentation rather than through the inter- 
mediary of a chart recording of the standing-wave pressure pattern. The second objec- 
tive was to increase the precision and accuracy of the data by incorporating a few simple 
and inexpensive techniques into the measurement procedure. A complete description of 
the experimental setup and test procedure is given in the following sections. 

Apparatus and Instrumentation 

Impedance tube.- A block diagram of the impedance tube apparatus and associated 
instrumentation is shown in figure 17. The apparatus consists of a main cylindrical tube 
section in which a test specimen is mounted at one end by means of a test fixture that 
provides an airtight seal around the specimen periphery. A variable cavity backing 
depth is provided by a solid movable piston in a cylindrical tube attached to the test fix- 
ture as shown in the sketch. The movable piston, 6.75 cm long, was machined from 
steel and was fitted with O-ring seals at each end to provide an airtight sliding contact 
with the machined inside surface of the cavity backing tube. Both the cavity backing tube 
and the main tube have an inside diameter of 5.72 cm and a wall thickness of 0.64 cm. 

The cavity backing tube and main tube have respective lengths of 60.95 cm and 83.83 cm. 

A sound source consisting of a 60-W electromagnetic driver was coupled to the main 
tube through an offset exponential horn as depicted in figure 17. A flexible coupling 
decoupled mechanical vibration of the driver from the tube walls. A further precaution 
was taken to reduce mechanical vibration of the tube walls by wrapping two layers of 
asphalt-based damping tape over the entire length of the tubes. 

Test specimen fix time. - The test specimen fixture was built to allow mounting of 
test specimens 6.22 cm in diameter and up to 1.27 cm thick. Specimens less than 
1.27 cm thick were bonded to an aluminum ring of sufficient length and thickness to pro- 
vide an axial clamping force around a peripheral strip of the specimen. An active diam- 
eter (i.e., diameter of specimen exposed to incident sound) of 5.72 cm remained. The 
outside surface of the aluminum mounting ring was coated with a thin layer of silicone 
release agent before installation to further insure an airtight seal between the incident 
and backii^ cavity side of the specimen. 

Acoustic pressure transducers and associated hardware.- Acoustic pressure signals 
could be monitored simultaneously at two of three locations by condenser microphones 
as shown in figure 17. A microphone 0.64 cm in diameter (channel 1) was used to 
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measvire the acoustic pressure level at the piston face. The velocity throt^h the speci- 
men could be calculated (thin specimens only) from this pressure level. A 1.27-cm 
microphone (channel 2) was coupled to a movable steel probe tube (O.D. = 0.33 cm; 

I.D. = 0.19 cm; Length = 122 cm) to obtain relative acoustic pressure level measurements 
at points on the axis of the main tube. Microphone 2 was also isolated mechanically from 
the probe tube and tube support hardware. Fixed positions of both the axial probe and 
piston microphone could be read to four significant figures (i.e., ±00.01 cm) with the use 
of a specially constructed vernier, which was rigidly attached to the probe tube. The 
axial probe could also be continuously traversed by means of a variable speed motorized 
gear box for recording standing -wave patterns with a graphic level recorder if desired. 

A third microphone 0.33 cm in diameter (channel 3) was mounted so that its axis coin- 
cided with a radius of the main tube. This microphone was positioned 0.69 cm from the 
incident side of the specimen; it could be traversed manually in the radial direction up to 
the tube axis. Also the separation distance between the radial microphone probe and the 
specimen surface coxold be changed by an amount equal to the specimen mounting ring 
length by reversing the specimen orientation. The geometry of this apparatus together 
with three channels of acoustic measurement capability permitted measurements of 
standing-wave ratios, standing-wave pressure minima positions, and acoustic pressure 
variations across the specimen face. 

Instrumentation. - A block diagram of the electronic instrumentation used in this 
experiment is also shown in figure 17. To insure reproducible measurement conditions 
at the test specimen, the acoustic pressure level was maintained at a constant level at 
the cavity back by monitoring the output of microphone 1. The signals from both the 
piston microphone and the axial probe microphone were filtered using a dual-channel 
10 -Hz -bandwidth tracking filter whose center frequency was automatically set by the 
oscillator frequency used to drive the acoustic source. A spectrum analyzer was used 
to monitor the overall spectral content of the signal at the piston face in order to detect 
malfunctions of equipment and/or nonlinear behavior of the acoustic source or test spec- 
imen. Generally, the harmonic content of the acoustic velocity at the specimen surface 
was maintained at about 20 dB below the fundamental acoustic velocity. The tracking 
filter outputs were read out visually on a dual-channel log-voltmeter which enabled 
acoustic pressure levels to be read to at least three significant figures over the entire 
test range. 


Test Specimens 

A summary of the geometrical and physical properties of the three test specimens 
used in this investigation are given in table I. The specimens were selected to cover the 
rai^e of acoustic impedance of most interest in aircraft duct liner applications. In 
addition, the fiber metal specimens B and C were selected to show the contaminating 
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effects, if any, of mechanical resonance on the intrinsic acoustic Impedance of the spec- 
imen, Specimen A was chosen as a control specimen in which the effects of mechanical 
vibration would be minimal and whose acoustic impedance could be computed theoretically. 

Specimen A.- Test specimen A was a honeycomb-like ceramic structure with 
a porosity of 59 percent (195 channels/cm^) and a thickness of 1.27 cm. The aver- 
age individual channel diameter was 0.062 cm. This specimen was characterized by 
a very high flexural stiffness and low bulk density. These mechanical properties 
insured that no mechanical resonance of the specimen would contaminate the acoustic 
impedance in the frequency range of interest. Also the acoustic impedance of this 
specimen could be predicted on the basis of the classical theory for sound propagation 
in a tube (ref. 8). The nominal flow resistance for this specimen was 2.52 cgs ray Is 

^1 rayl = 1 ■ = 10 Ns/m3^ . 

Specimen B.- Test specimen B consisted of fiber metal fabricated from fibers of 
AISI type 347 stainless steel in a random arrangement. The fibers were compressed 
to a thickness of 0.09 cm. The mechanical properties of the specimen were such that 
the first mechanical resonance occurred within the test frequency range. The nomi- 
nal flow resistance for this specimen was 57 cgs rayls. 

Specimen C.- Test specimen C was also fabricated from fibers of AISI type 347 
stainless steel compressed to a thickness of 0.35 cm. The nominal dc flow resistance of 
this specimen was 420 cgs rayls. The flexural stiffness and density were such that two 
mechanical resonances occurred within the range of frequencies investigated. 

Tube Checkout Procedure for Systematic Errors in Reflection 
Coefficients and Pressure Minima Positions 

In addition to tube-wall absorption effects, systematic errors in the reflection coef- 
ficient and pressure minima positions can result from 

(a) Excitation of higher order acoustic modes in the impedance tube caused by tube 

axis curvature at the probe entrance (see fig. 17) or other tube irregularities 

(b) Departure of the backing cavity piston from assumed perfect reflectivity due to 

possible air leaks, microphone installation, vibration, thermal bovindary 

layer, etc. 

(c) Flanking transmission of sound through axial probe walls, acoustic excitation of 

axial probe, or other mechanical disturbances 

(d) Excessively high ambient noise floor of acoustic medium in the impedance tube 

(e) Electrical noise floor of measuring equipment 
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Unlike tube-wall absorption effects, possible systematic errors due to the above sources 
must be evaluated experimentally. The overall effect of the factors just listed can be 
estimated by observing the wave field in the tube in the absence of a test specimen; the 
backing cavity piston should be placed at the test specimen location (see fig. 17). This 
procedure is accomplished by observing amplitude and phase variations in the radial 
direction and by measuring reflection coefficients and pressure minima positions. 

Radial amplitude and phase variations.- The radial microphone probe was manually 
traversed across about 75 percent of the tube diameter at each test frequency. The 
relative amplitude and phase changes between the axial microphone probe 2 and radial 
microphone probe 3 were observed with the axial probe fixed. Over the 0.5- to 3.0-kHz 
frequency range, the relative amplitude and phase changes were undetectable. Over the 
3.0- to 3.5-kHz frequency range, the relative amplitude and phase changes were approxi- 
mately 0.2 dB and 2°, respectively. Above 3.6 kHz the relative amplitude and phase 
changes were clearly indicative of the existence of nonplanar waves. The theoretical 
cuton frequency of the first axisymmetric higher order mode was approximately 3.6 kHz. 

Reflection coefficient behavior with no test specimen.- An additional check for 
anomalous behavior of the wave field in the impedance tube apparatus as well as the 
establishment of the signal noise floor was accomplished by measuring the standing- 
wave ratios and corresponding pressure minima positions with no test specimen present 
and with the piston face coincident with the specimen reference plane. This procedure 
also established the magnitude of possible systematic errors associated with the pressure 
minima positions and reflection coefficients. Figure 18 shows best estimates of the 
reflection coefficients and the corresponding best estimates of the nondimensional first 
pressure minima positions for frequencies of 0.5 to 3.5 kHz (at intervals of 0.1 kHz) with 
the SPL at the piston face held constant at 120 dB. The formula of equation (1) was used 
to calculate the reflection coefficients from the measured standing -wave ratios. These 
coefficients are plotted in figure 18 along the right-hand ordinate. These data show no 
anomalous behavior in the tube over the indicated frequency range. Although the reflec- 
tion coefficients vary with frequency, they do not drop below 0.9950, and they remain 
above 0.9980 for most frequencies. These fluctuations are associated with the large 
standing -wave ratios in the tube as indicated by the dashed horizontal lines. Note that 
an Sj^ of 60 dB implies an SPL of 60 dB at a pressure minimum if the level at the piston 
face is 120 dB. These low SPL regions occupy very thin cross sections of the tube. 
Hence, any of the items listed under (c), (d), or (e) in the discussion of tube checkout pro- 
cedure can cause an increase in the signal noise floor with an accompanying apparent 
reduction in the standing-wave ratio. This conjecture was supported by the fact that 
fluctuations by as much as 20 dB occur for the higher values of S^^. These measure- 
ments indicate that the piston behavior approaches that of a perfect reflector over the 
frequency range 0.5 to 3.5 kHz. 
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Pressure minima behavior with no test specimen.- For a perfectly rigid reflector 
the first pressure minima for all wavelengths should cluster around Xj = 0.25 provided 
random measurement errors only are present. However, the mean trend of the first 
pressure minima in figure 18 suggest a constant systematic error which would result in 
a displacement of the pressure minima positions 0.09 cm too far from the test specimen 
surface. As discussed previously in the subsection entitled "Systematic Errors," this 
discrepancy can be caused by a combination of scale calibration error and finite probe 
size effects, both of which are frequency independent. The systematic error in pressure 
minima positions was removed by incorporating a correction into the computer program 
used for the computations discussed previously in the "Data Processing" subsection. 

In summary, the data of figure 18 demonstrate that the movable piston arrange- 
ment with the flush mounted microphone approaches a perfectly rigid reflector in its 
acoustical behavior. However, corrections for systematic errors in measurements for 
the pressure minima positions were necessary. 

Tube calibratio n for random errors in reflection coefficients and pressure minima 
positions. - Random measurement errors occurring in the reflection coefficients and 
pressure minima positions were treated as properties of the impedance tube apparatus 
and its operating environment. As previously indicated, such treatment is a plausible 
procedure for moderate test specimen impedances. Hence, statistical data relevant to 
the random measurement errors can be combined as different specimens are tested. 
When sufficient statistical data have been collected to construct probability distributions, 
these distributions can be considered properties of this particular experimental arrange- 
ment just as was the case with systematic errors. 

At least three pressure minima positions and corresponding standing -wave ratios 
were measimed at each test frequency. Data collection efficiency was improved by 
recording pressure minima positions and SPL readings directly on data sheets from 
which computer inputs in the form of punched cards could be prepared. This procedure 
eliminated the need for time-consuming chart recordings of the entire standing-wave 
pattern. Accuracy in measurement of pressure minimum position was improved by tak- 
ing two scale readings symmetrically spaced on either side of a given pressure mini- 
mum. A pressure minimum position was then obtained by taking the average of these 
readings. 


RESULTS AND DISCUSSION 
Impedance Measurement Results 

Figures 19 to 26 present the experimental and theoretical results of this investiga- 
tion for the three test specimens designated A, B, and C. Resistance and reactance 
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ratios together with their calculated fractional standard deviations are presented for each 
of the three test specimens as a function of frequency. In addition, the effect of system- 
atic changes in the pressure minima positions is shown in these figures. For most of 
the tests, the backing cavity depth was adjusted to one-quarter wavelength for each test 
frequency, and the SPL at the cavity back was held constant at either 90 or 114 dB. The 
results of one test are presented in which the cavity depth was held constant. These 
results illustrate the possible systematic errors that may arise due to cavity antireso- 
nance. Finally, for specimens B and C, the effect of mechanical vibration of the test 
specimen on the intrinsic acoustic impedance of the specimen is illustrated. 

Specimen A 

Effect of sys_tematic measurement errors.- The effects of systematic errors in 
pressure minima positions and reflection coefficients were investigated by perturbing the 
measured quantities. Perturbations of approximately 5 times the respective standard 
deviations as obtained from the data of figures 3 and 5 were used. The effect of a pertur- 
bation of ±0.05 cm in Xj ^i.e., ^ is shown in figures 19 and 20 by the short-dashed 

curves that diverge about the measured impedance data as the frequency increases. The 
±0.05-cm perturbation was applied to the test data represented by the circular symbols. 

A similar calculation to evaluate the systematic errors due to a ±0.005 perturbation 
^i.e., ±5 ctj^^^ in the reflection coefficient indicated a relatively insignificant effect com- 
pared with systematic errors in x^^ and therefore were not plotted. These results illus- 
trate that the distance to the first pressure minimum expressed in wavelengths ^i.e., Xj^ 
is critical in the measurement of impedance values produced by specimen A. At 3.5 kHz, 
for instance, the resistance and reactance ratios can be increased by approximately 10 
and 6 percent, respectively, by a 1.5-percent increase in Xj (i.e., 0.05-cm increase 
in . Note, however, that the three sets of data shown in figures 19 and 20 do not show 
a systematic increase in the scatter with increasing frequency. Systematic error in the 
measurements does not appear to be a significant problem over the time span and test 
conditions for which these measurements were obtained. 

Effect of random measurement errors.- The effect of random measurement errors 

in pressure minima positions x^^ and reflection coefficients as specified by 

a = 0.0096 cm and Op = 0.001, respectively, are shown in figures 19 and 20 by the 
i 1 

long-dashed curves. These curves represent the fractional standard deviations in the 
resistance and reactance ratios, respectively, as a function of frequency. The curves 
were calculated using the approximate method as given in equations (18) to (21). In all 
cases, two sets of mean values for resistance and reactance were used in equations (18) 
to (21), corresponding to the two short-dashed curves representing the effect of a 
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±0.05-cm perturbation on Xj^. This was done to reveal the effect of mean value changes 
in resistance and reactance ratios on their respective fractional standard deviations. 

For specimen A, figure 19 shows that no significant chaise in the resistance ratio frac- 
tional standard deviation occurs for the changes in the mean values of the resist- 

ance ratio 6 associated with the short-dashed curves. Also, for the reactance ratio 
standard deviation a /\ the mean value changes corresponding to the short-dashed 
curves of figure 20 w^e found not to have a significant effect; therefore, only one curve 
was plotted for Examination of the magnitudes and trends of 

plotted against frequency in figures 19 and 20, respectively, indicates a generally 
increasing effect of random errors with increasing frequency. However, and 

cr exceed approximately 2 and 1.3 percent, respectively. Specifically, at 
3.4 kHz, there is a 67-percent probability that the measured resistance ratio lies within 
approximately 0.278 ± 0.006 pc units. The range of variation corresponds closely to the 
amount of scatter observed in the three measurements at this frequency. For the 
reactance ratio the worst scatter appears to occur at 3.3 kHz where ^y^X is about 
1.1 percent. Therefore, there is 67-percent probability that the reactance ratio lies 
within approximately 1.52 ± 0.02 pc units. In this case, the observed data scatter is 
slightly greater than 2a^. These observations suggest that the contribution of random 
errors in pressure minima positions and reflection factors to the resultant scatter in the 
deduced impedance components is adequately described by the approximate calculation 
procedure given by equations (18) and (21) if the limits of a and discussed in 

Xi K. 

figures 15 and 16 are not exceeded. 

Comparison of predicted and measured impedance.- Of the three test specimens 
only the surface impedance of specimen A could be predicted from existing classical 
theory relating to sound propagation throi^h small tubes (ref. 8). An elementary appli- 
cation of the results of reference 8, together with the following assumptions, provides a 
predictive model for the surface impedance of specimen A backed by a cavity with an 
impedance 

(1) The propagation constant T for the entire length of a channel is equal to that 
for an infinitely long channel of the same radius. 

(2) All channels have uniform circular cross sections of equal radii. 

(3) The solid structure containing the channels is perfectly rigid. 

(4) The ratio of channel diameter to channel separation is sufficiently near unity to 
eliminate the effects of the classical Rayleigh end correction (ref. 9). 

With these assumptions, the surface impedance of a high porosity structure with parallel 
channels is given by 
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(22) 



PZ^j + tanh n 
P \l + tanh nj 


where Z^ is the impedance of the backing cavity, Z^ is the characteristic impedance 
of a single channel in the specimen, P is the specimen porosity, P is the propagation 
constant for waves propagating through a channel, and I is the channel length. If the 
channel backing cavity depth is set equal to a quarter-wavelength, as was the case in the 
investigation, then 7i-^ = 0 and equation (22) reduces to 


Z = — tanh n (23) 

P 

The propagation constant is given in reference 8 as 



(24a) 


where s is the shear wave number. The quantity n is a complex polytropic constant 
describing the thermodynamic process in the channel and is given by reference 8 as 


n = 



(24b) 


where /3 is the square root of the Prandtl number. The characteristic impedance is 
given by 



(24c) 


2 

The shear wave number s and the Prandtl number ^ are nondimensional quantities 
and are given, respectively, by 
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(24d) 


s = r 



/3 = 



(24e) 


The geometrical properties of specimen A needed in equations (23) and (24) are given in 
table I. In addition, the appropriate physical properties of the air in the channels of 
specimen A were talcen to be 

T = 25° C 

p = 0.0018 g-cm“^ 

ji^ = 0.000181 g-cm"l-sec"l 

c = 34 645 cm/sec 


/3 = 0.88 
y = 1.402 

The continuous curves of figures 19 and 20 show the predicted values of resistance 
and reactance ratios, respectively. The theory accurately predicts the measured resist- 
ance ratio at the lower frequencies but overpredicts the measured values by as much as 
24 percent at the higher frequencies. For the reactance ratio, the theory overpredicts 
the measured values by about 14 percent at the lower frequencies; the overprediction 
increases to about 24 percent at the higher frequencies. A sensitivity study of equation (23) 
indicated that the predicted resistance and reactance ratios were much more sensitive to 
changes in the channel geometry than to changes in the gas temperature. These results 
suggest that the assumptions regarding end effects and channel geometry may be inad- 
missible for structures similar to specimen A; however, derivation of a more realistic 
model was beyond the scope of this investigation. 

Specimen B 

Specimen B consisted of felted stainless -steel fibers compressed into a thin (rela- 
tive to the acoustic wavelength) sheet. The physical properties of the test specimen are 
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given in table I. This material has been a contender for the resistive element in air- 
craft duct liner applications. This particular specimen was chosen for this investigation 
because the resistance was such that near optimum performance of the standing -wave 
method for measuring impedance could be expected (i.e., the standing -wave ratios ranged 
between 3 and 20 dB). 

Test results (quarter -wavelength backing cavity depth) .- Figures 21 and 22 show the 
measured resistance and reactance ratios, respectively, for five sets of impedance data 
for test specimen B. These data were collected over a time span of approximately 
9 months. Also, the specimen was removed and reinstalled between tests. Following the 
same format used in figures 19 and 20 for specimen A, a ±0.05-cm perturbation of the 
pressure minima positions for a typical test run (the impedance data collected on 
April 19, 1976, in this case) are shown by the short-dashed curves. Also, the long-dashed 
curves represent the behavior of the fractional standard deviations according to the mean 
values taken from the short- dashed curves. 

Except for the sharp decrease in the resistance and reactance ratios in the vicinity 
of 1.4 kHz, the measured impedance values for all five sets of data exhibit uniform scatter 
across the entire 0.5- to 3.5-kHz frequency range. The dramatic decrease of the resist- 
ance and reactance ratios in the vicinity of 1.4 kHz results from the fundamental mechan- 
ical resonance frequency of the test specimen. This result was confirmed by calculating 
the resonance frequency based on estimates of the mechanical stiffness and surface den- 
sity of the test specimen. (It should be noted that a cavity backing depth of a quarter- 
wavelength gives rise to the maximum possible acoustic pressure across the specimen 
surface.) These data provide an opportunity to better understand the relative contribu- 
tions of systematic and random measurement errors. 

As illustrated by the resistance ratio data in figure 21, the small percentage change 
(0.3 to 1.1 percent) in the fractional standard deviation over the 0.5- to 3.5-kHz frequency 
rar^e does not adequately account for the generally uniform scatter in the data (outside 
the resonance region). On the other hand, systematic errors (associated with pressure 
minima measurements) between the different sets of measurements should produce 
diverging resistance values with increasing frequency as illustrated by the short-dashed 
curves. (However, no significant divergence was found for the reactance ratio when the 
measured positions were perturbed by ±0.05 cm.) These observations suggest that neither 
systematic errors intrinsic to the impedance tube nor random measurement errors are 
significant contributors to the observed scatter. The remaining likely source of data 
scatter was systematic error associated with the test specimen movinting conditions. 

This possibility can be supported by observing the detailed behavior with frequency of 
individual tests, especially in the vicinity of the mechanical resonance. In particular, 
note that the resistance ratio changed from 1.18 at 1.45 kHz for the test conducted on 
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April 19, 1976, to 0.95 at 1.45 kHz for April 21, 1976, or a change of 19 percent over a 
2 -day time period. Similar changes can be observed at other frequencies in the neighbor- 
hood of the mechanical resonance. At frequencies removed from the region of resonance, 
the variation of the measured resistance ratio among the different tests is approximately 
a uniform 4 percent. 

The changes in the measured impedance values shown in figures 21 and 22 appear 
to be random in a manner independent of frequency. These observations suggest that the 
mounting conditions (i.e., conditions at the clamped boundary of the specimen) are chiefly 
responsible for the observed differences between the repeated impedance measurements. 

It should be noted that the axial clamping force around the specimen boundary was uni- 
formly applied by means of machine bolts (1/4 in. diam.; 28 threads/in.) on which torque 
was applied until the specimen test fixture flange was seated against the impedance tube 
flange. This procedure was designed to allow the same amount of specimen boimdary 
compression for each test run. The torque on the bolts was not necessarily the same 
from test to test. 

Fixed cavity backing depth .- The design of practical duct liners for noise reduction 
applications may involve several layers of resistive face-sheet material together with 
various cavity depths to provide desired broadband impedance behavior. It is important 
to verify experimentally the predicted impedance of such structures. Therefore, the 
limitations of the particular impedance measurement method used to verify theoretical 
predictions should be understood. The transmission line method appears to have a funda- 
mental limitation which will be illustrated with the results of a simple test. Specimen B 
was provided with a fixed cavity depth corresponding to a quarter -wave length at 1.0 kHz. 
The impedance of the cayity-specimen combination was then measured over the 0.5- to 
3.5-kHz frequency range with the sound pressure level at the cavity back held constant 
at 90 dB. 

Figures 23 and 24 show the results of the fixed cavity depth impedance measure- 
ments for specimen B. When the backing cavity is a pure reactance, the measured 
resistance should theoretically be that caused by the face sheet alone as was the case 
for the variable quarter -wave length cavity. However, when the acoustic half-wavelength 
is equal to the cavity depth, the acoustic particle velocity tends toward 0 at the specimen 
face. Physical intuition would suggest that the impedance would become very large since 
the acoustic pressure is held constant. The result illustrated in figure 23 for the resist- 
ance ratio plotted against frequency indicates that the intuitively suspected large imped- 
ance is achieved by the resistance becoming very large. Note that the theoretical cavity 
reactance (solid curve in fig. 24) becomes indeterminant at 2.0 kHz. The measured 
reactance closely follows the cavity reactance at low frequencies and becomes slightly 
greater than the cavity reactance with increasing frequency except in the vicinity of 
2.0 kHz. This general deviation of the measured reactance toward more positive values 
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is due to the contribution of the positive reactance of specimen B to the cavity reactance. 
Since the measured reactance should theoretically never be less than the cavity reactance, 
the measured reactance at frequencies 1.9 kHz and 2.0 kHz in figure 24 must be regarded 
as anomalous. The possible effects of random measurement errors are indicated by the 
dashed curves in figures 23 and 24. Note that the fractional standard deviation rises 
dramatically near 2.0 kHz. At the frequency of 2.0 kHz where the maximum possible 
effect of random measurement errors occurs, the 2 ct^ and 2og ranges have been indicated 
on the measured impedance components in the figures. Although there is significant uncer- 
tainty in the impedance components due to random measurement error, the measured sys- 
tematic increases in resistance are not fully explained by this uncertainty. 

The results of the fixed cavity depth impedance measurements depicted in fig- 
ures 23 and 24 suggest that the method does not produce valid results in the vicinity of 
those frequencies that produce multiple half-wavelengths in the backing cavity. The 
anomalous behavior of the measured impedance may be due to some as yet undetermined 
systematic error in the experimental setup. 

Specimen C 

Specimen C consisted of felted stainless steel fibers compressed into a thin (rela- 
tive to the acoustic wavelength) sheet. The physical properties of the specimen are given 
in table I. Specimen C was chosen for this investigation because of its relatively high 
flow resistance of 420 cgs ray Is. A flow resistance of 420 cgs rayls is as high as is 
likely to be encountered in any practical duct liner application. A material with such 
high resistance produces pressure minima positions in the vicinity of the critical quarter- 
wavelength point relative to the specimen face. These positions should reveal the relative 
importance of random and systematic errors for highly resistive specimens in view of 
the behavior of the error functions Fj and F 2 (see figs. 9 and 10). Therefore, spec- 
imen C should serve to define the importance of random and systematic errors in the 
measured impedance for highly resistive specimens. 

Figures 25 and 26 show the measured resistance and reactance ratios, respectively, 
for five sets of data collected over a time span of approximately 1 month. All the data 
were collected with a variable cavity backing depth of a quarter-wavelength. The speci- 
men was removed and reinserted between tests. Following the same format used for 
specimens A and B, the effects of systematic errors of ±0.05 cm in the data for pressure 
minima positions collected on March 8, 1976, is shown by the short-dashed curves. The 
long-dashed curves represent the behavior of the fractional standard deviations using the 
mean values taken from the short-dashed curves. 

A number of interesting features are evident from figures 25 and 26. First, two 
mechanical resonances contaminate the intrinsic acoustic impedance of the test specimen. 
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These resonances occur at approximately 0.67 and 2.15 kHz and correspond to the funda- 
mental and first harmonic of the specimen acting as a clamped plate. Calculated reso- 
nance frequencies based on estimates of the specimen stiffness and area density indicated 
that the second resonance frequency corresponded to a vibrational mode with one circum- 
ferential node line. Second, in contrast to the situation for specimens A and B, the frac- 
tional standard deviations were significantly sensitive to the mean values of the impedance 
ratios taken from the short-dashed curves. For the case of the reactance, the 0.05-cm 
perturbation shifts the reactance to the vicinity of 0 over the frequency ranges from 1.0 j 

to 2.0 kHz and from 2.5 to 3.4 kHz. The resulting effect on the reactance fractional 
standard deviation is meaningless. The corresponding long-dashed curve indicates that 
the fractional standard deviation exceeds 120 percent in these regions. The reactance 
data appear to support these analytical results since the reactance ratios exhibit more 
scatter with increasing frequency than do the resistance ratios. Third, the reactance 
ratios are critically sensitive to systematic error in the locations of pressure minima 
positions at the higher frequencies (see the short-dashed curves). In particular, careful 
examination of data for pressure minima positions for March 2, 1976, reveals that xj 
experienced a 0.0023 increase (0.082 cm) between 2.7 and 2.8 kHz. This shift caused the 
sudden increase in the reactance ratio measured on that date. This sudden increase 
deviated from the pressure minima position trends for the other tests in which the pres- 
sure minima positions were moving toward the test specimen surface. A similar, 
apparently inconsistent, variation in pressure minimum position caused the sudden drop 
in the reactance ratio at 1.8 kHz for the data taken on March 3, 1976. Furthermore, 
corresponding disturbances in the resistance ratio trends are less severe. A careful 
review of the tube performance data did not reveal any significant inconsistencies. 

Therefore, these anomalies in the impedance must have resulted either from small dis- 
turbances in the probe position calibration or from contaminating influences of the spec- 
imen’s mechanical impedance on the specimen's intrinsic acoustic impedance. 

CONCLUSIONS 

An investigation by the transmission line method of the effect of random and system- 
atic errors associated with the measurement of normal incidence acoustic impedance in 
a zero-mean-flow environment indicated the following conclusions: 

1. The effect of random measurement errors on impedance as measured by the 
transmission line method can be conveniently evaluated from estimates of the variances 

and O'r ^ where Xj is the distance to the first pressure minimum and Rq is 

the test specimen reflection coefficient. 
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2. Comparison of an exact and an approximate calculation procedure for obtaining 

2 2 

the variances of the resistance and reactance ratios Cn and ct, indicated that com- 

p X 

putation times could be reduced by a factor of 100 without significant reduction in accu- 
racy provided that standard deviations a and CTd do not exceed 0.045 cm and 0.006, 

Xj n. 

respectively. The approximate formulas are given. 

3. Random measurement errors in pressime minima positions and reflection coef- 
ficients for the experimental setup used in this investigation were found to approximate 
a normal probability distribution with standard deviations of 0.0096 cm and 0.001, 
respectively. 

4. Calculated fractional standard deviations in the measured impedance components 
were generally fovind to increase with frequency; however, sensitivity to mean values of 
the measured quantities was found to vary markedly. 

5. Conclusions 3 and 4 led to the belief that the uniform scatter in the impedance 
components observed to occur with increasing frequency for repeated measurements 
resulted from variations in test specimen mounting conditions. 

6. Impedance measurements are, in general, more sensitive to random or system- 
atic errors in the pressure minimum position than to random or systematic errors in the 
reflection coefficients. 

7. Reliable resistance measurements can be obtained with constant cavity depth 
testing except in the vicinity of frequencies that produce multiples of a half-wavelength 
in the cavity. 

These conclusions suggest that the transmission line method for measuring acoustic 
impedance is attractive for its theoretical and operational simplicity. Furthermore, the 
method allows the effects of random measurement errors (precision) to be quantified with 
ease. For the range of impedance values of interest in aircraft duct liner applications, 
pressure minimum position and standing -wave ratio measurements can easily be refined 
to limit impedance data scatter intrinsic to the tube to insignificant amounts (standard 
deviations not greater than 2 percent of the mean below 3.5 kHz). Careful consideration 
of the impedance measurements and error propagation analysis indicate that mechanical 
bovmdary conditions imposed on the specimen constitute the major source of inaccuracy 
in the measurement of isolated face-sheet impedance values. These boundary conditions 
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not only cause severe contamination of the intrinsic acoustic impedance of the test spec- 
imen in the neighborhood of mechanical resonances but may also cause most of the vari- 
ation of the impedance components between identical tests. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 30, 1977 
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TABLE I. - GEOMETRICAL AND PHYSICAL PROPERTIES OF TEST SPECIMENS 


Specimen 

designation 

Description 

Thickness, 

cm 

Active 

diameter, 

cm 

Bulk 

density, 

g/cm3 

Porosity, 

percent 

Fiber 

diameter, 

cm 

Channel 

diameter, 

cm 

Flow resistance, 
cgs rayls 
(at 1 cm/sec) 

A 

Ceramic 

honeycomb 

1.27 

5.71 

0.62 

59.0 - 


0.062 

2.52 

B 

Fiber metal 

0.09 

5.71 

5.40 

29.9 

0.0112 

— 

57.0 

C 

Fiber metal 

0.35 

5.71 

1.81 

76.5 

0.0004 

— 

420.0 




r 


4 :^ 

O 



X, cm 


Figure 1.- Impedance tube and typical standing -wave pattern. 





X, cm 

Figure 2.- Sketch showing correspondence between standing -wave 
ratios and reflection coefficients with linear regres- 
sion curve to obtain Rq. 
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e to pressure minima 












Figure 9.- Error propagation function Fj plotted against distance to first pressure 
minimum in wavelengths with reciprocal standing-wave ratio rg as parameter. 
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Error propagation function, 





Fractional standard deviation, > percent 



(a) Resistance ratio. 

Figure 11.- Fractional standard deviation in impedance ratio components plotted against reflection 
coefficient with distance to first pressure minimum as parameter for frequency of 0.5 kHz. 
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(b) Reactance ratio. 
Figure 11.- Concluded. 
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Reflection coefficient, Rq 
(a) Resistance ratio. 

Figure 12.- Fractional standard deviation in impedance ratio components plotted against reflection 
coefficient with distance to first pressure minimum as parameter for frequency of 1.0 kHz. 
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(a) Resistance ratio. 

Figure 13.- Fractional standard deviation in impedance ratio components plotted against 
reflection coefficient with distance to first pressure minimum as parameter for fre- 
quency of 2.0 kHz. 
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(b) Reactance 
Figure 13.- Co: 
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Reflection coefficient, Rq 


(a) Resistance ratio. 

Figure 14.- Fractional standard deviation in impedance ratio components plotted against 
reflection coefficient with distance to first pressure minimum as parameter for fre- 
quency of 3.0 kHz. 
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Figure 15.- Fractional standard deviation in impedance ratio components plotted against 
standard deviation in reflection coefficients with distance to first pressure tnlnimnm 
as parameter. 
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(a) Resistance ratio. 

Figure 16,- Fractional standard deviation in impedance ratio components plotted against 
standard deviation in distances to pressure minima with distance to first pressure 
minimum as parameter. 
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Figure 17.- Block diagram of experimental setup. 
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